hep-lat/96l'Jum 



LSUHE No. 249-1996 
HD-THEP-96-58 
ILL-(TH)-96-17 
OUTP-96-75P 
December, 1996 



the Phase Diagram of the SU(2) Adjoint Higgs Model in 2+1 Dimensions 



A. Hart 1 , 0. Philipsen 2 , J. D. Stack 3 and M. Teper 4 

Department of Physics and Astronomy, Louisiana State University, 
Baton Rouge. LA 70803 U.S.A. 

2 Institut fur Theoretische Physik, Universitat Heidelberg, 
Philosophenweg 16, D-69120 Heidelberg, Germany. 

^Department of Physics, 
University of Illinois at Urbana- Champaign, 
1110 W. Green Street, Urbana. IL 61801 U.S.A. 

4 Theoretical Physics, University of Oxford, 
1 Keble Road, Oxford. 0X1 3NP U.K. 

(to be published in Physics Letters) 



Abstract 



The phase diagram is investigated for SU (2) lattice gauge theory in d = 3, coupled to 
adjoint scalars. For small values of the quartic scalar coupling, A, the transition separating 
Higgs and confinement phases is found to be first-order, in agreement with earlier work 
by Nadkarni. The surface of second-order transitions conjectured by Nadkarni, however, is 
shown instead to correspond to crossover behaviour. This conclusion is based on a finite 
size analysis of the scalar mass and susceptibility. The nature of the phase transition at 
the termination of first-order behaviour is investigated and we find evidence for a critical 
point at which the scalar mass vanishes. The photon mass and confining string tension are 
measured and are found to be negligibly small in the Higgs phase. This is correlated with 
the very small density of magnetic monopoles in the Higgs phase. The string tension and 
photon mass rise rapidly as the crossover is traversed towards the symmetric phase. 



In this paper we report on our investigation of the phase diagram for d = 3 SU(2) 
lattice gauge theory coupled to adjoint scalars. This theory is of interest for at least two 
reasons. Through dimensional reduction, it is related to d = 4 pure SU (2) lattice gauge 
theory at high temperature. The model also contains in its Higgs phase the famous 't 
Hooft-Polyakov monopoles. The Higgs phase is continuously connected to the symmet- 
ric, confining phase. Thus the model may prove to be a useful laboratory for exploring 
confinement by monopoles. We will not further discuss monopole confinement here, but 
concentrate on the phase diagram and some general properties of the mass spectrum. 

The continuum Euclidean action in the notation of Nadkarni is given by 

S = Jd 3 x j^Tr(F^) + TriD^D^) + m^Tr(#) + |(Tr(#)) 2 | , (1) 

where F^ v = d^A u - d v A VL + i[A^A v ], D^<p = d^<j) + iiA^tfi]; and F^A^ and <p are all 
traceless 2x2 Hermitian matrices {(f) = <p a <j a /2, etc). The lattice form of the action is 

S = P E fl-iTr[/ Ml/ (x))+2^Tr($(x)$(x)) 



X,[t>l/ s ' X 

- 2/€^Tr($(x)^(x)$(x + /ia)f/ t (x)) + A^(2Tr($(x)$(x)) - l) 2 (2) 

X,fl X 

Here U^(x) = exp(iagA^(x)) is the usual link variable, and U^x) is the corresponding 



plaquette. The continuum field 4>{x) is related to the lattice field by <f>(x) = y/t/a$(x), 
where a is the lattice spacing, and the lattice parameters (3 and A are defined by 

_4_ ™£ 
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The hopping parameter k is determined by an equation relating it to the continuum mass 
parameter m\. The tree level form of this equation is 
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Dividing by a factor g 4 , and using Eq.(D), this can be rewritten as 



Due to the superrenormalisability of the adjoint Higgs theory in d — 3, A and g 2 do not 
require ultraviolet renormalisation. Suppose for the moment that this were true for ttIq 
as well. Eq.(^j) then specifies how to adjust k to correspond to given fixed ratios of the 
continuum parameters m^j 'g 4 and p/g 2 . The lattice spacing is set by a = 4/g 2 f3. 

The mass parameter Wq does of course require ultraviolet renormalisation, but only up 
to two loops. The known results in the MS scheme [[§ and in lattice perturbation theory 
H allow a generalisation of Eq. (||) to be derived, 

m 2 {g 2 ) 

9 4 ~ 8 [k Pg 2 ) 
4vr 1 8g 2 



1 

+ 



167T 2 



(5) 



where E = 3.17591. The successive lines of this equation correspond to tree level, one loop 
and two loop terms, respectively. The dimensionful parameter in the MS scheme has been 
chosen to be g 2 . Eq.(H) specifies how to approach the continuum limit (/3 — > oo, A — > 0, k — > 
1/3), with m^{g 2 )/g i and \xj ' g 2 held fixed. It becomes exact in the limit of vanishing lattice 
spacing. 

Our analysis of the phase diagram is based on the early work of Nadkarni [I]], who first 
investigated this model using a combination of numerical and analytic techniques. Following 



Nadkarni, we define the length p of the scalar field by p = y2Tr($$). Qualitatively the 
Higgs phase and symmetric phases can be distinguished by the value of (p) , it being larger in 
the Higgs phase. The familiar language of symmetry breaking is useful in the Higgs phase, 
whereas the symmetric phase is better described as a confined hadronic phase. As was 
emphasized by Nadkarni, however, the two phases are in fact continuously connected, and 
there is really no fundamental distinction between them. When an actual phase transition 
separates the two phases, it is described by a surface « c (/3, A), the Higgs (symmetric) phase 
lying above (below) the surface. The statement that the two phases are continuously 
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connected implies that this surface does not cover the entire (3 — A plane, so that any two 
points in the parameter space can be connected by a path which detours around it. While 
it is not feasible to map out the entire surface in a strictly numerical analysis, one of our 
principal results will be that the surface of phase transitions covers a much smaller region 
of the (3 — A plane than conjectured by Nadkarni. As discussed in more detail below, the 
vast region of the surface specified by Nadkarni as second-order corresponds rather to a 
crossover region which is without critical fluctuations or massless particles. 

Our methods are similar to those used by two of the present authors in a previous study 
in d = 3 with SU(2) gauge fields coupled to fundamental rather than adjoint scalars Pj. 
The link update used in M was modified by inserting a Metropolis step to take account 
of the quadratic dependence of the scalar term in Eq.(|2]) on the link variables. The Bunk 
algorithm for the scalar update was used here also. With adjoint scalars, however, 
the approximate determination of Bunk's parameter a suggested in || and used in Q 
occasionally resulted in a low acceptance rate for the scalar update (particularly at very 
low (3). Solving the full cubic equation (Eq.(6) of f5|) for a gave acceptance rates of over 
80% throughout the phase diagram. 

The results presented here were obtained on lattices of size L 3 with L = 8, 12, 16, 20, 
24 and 32 at (3 = 6.0 and 9.0. A typical run consisted of 500 - 2000 sweeps. The coupling 
(3 = 6.0 is already in the region of approximate scaling for d = 3 pure SU(2) lattice gauge 
theory ||. In the case with fundamental scalars 0], the mass spectrum at the nearby value 
of 13 = 7.0 has been calculated and has been found to scale very well to higher values of 
(3. In the present case of adjoint scalars, we shall see later that there is good scaling of 
mass ratios between (3 = 6.0 and 9.0. The procedure is, given k and A at say, (3 = 6.0, to 
use Eqs.@ and ([5]) to determine the ratios of continuum parameters m^(g 2 )/ ' g A and /i/g 2 . 
Then holding these ratios fixed, Eqs.(|3|) and @ are solved again to find k and A at the 
new value of (3, e.g. (3 = 9.0. Close enough to the continuum limit, physical mass ratios 
should be constant when moving along such a 'line of constant physics'. 



3 




Figure 1: Probability distributions of p indicating a first- order phase transition at (3 = 6.0. 
(a) k = 0.370473, A = 0.001, L = 8. (b) k = 0.38887, A = 0.007, L = 20. 

The most straightforward region of the phase diagram to analyse is where Higgs and 
symmetric phases are separated by a first-order transition. In agreement with Nadkarni, we 
find clear first-order behaviour for small A, which for (3 = 6.0 means A <C 0.01. In Fig. [T|(a) 
we show the histogram of p at (3 = 6.0, A = 0.001 on an 8 3 lattice, p is the average over 
the lattice volume, V, of p(x). Two peaks corresponding to symmetric and Higgs phases 
are clearly seen. Tunnelling between the two states occurred once in this run. When we 
increased the lattice size the peaks became so distinct that the to-and-fro tunnelling became 
unobservably rare. The first-order nature of the transition can also be seen from a plot of 
(p) versus K. There is a jump from smaller to larger values of (p) which becomes rapidly 
sharper as the lattice volume is increased. The first-order behaviour observed at A = 0.001 
persists but weakens as A increases. The two peaks in p become less widely separated and 
tunnelling between the two states becomes more frequent, e.g. Fig. |l|(b) for A = 0.007. 
Here the barrier is weak and we have frequent tunnelling, which persists even on fairly large 
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lattices. Clearly at this value of A the first-order transition has become very weak and, as 
we shall see below, for slightly larger A one loses all evidence for first-order behaviour. We 
see exactly the same behaviour in the corresponding calculations at (3 — 9.0. 
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Figure 2: Finite size scaling of scalar susceptibility and mass indicating crossover be- 
haviour, (a) (3 = 6.0, A = 0.01. (b) (3 = 9.0, A = 0.01. (c) (3 = 6.0, A = 0.01. 

Let us defer for the moment the question of what happens at the precise value of A 
(for fixed (3) at which first-order behaviour terminates. It was conjectured by Nadkarni 
that beyond the region of first-order behaviour, the surface k c (/3, A) continues (all the way 
to A = oo at (3 = 6.0 and 9.0), describing a second-order transition. We instead find no 
actual phase transition, only a smooth crossover behaviour. This conclusion is reached 
using a finite size analysis of the scalar mass and related susceptibility. Our scalar mass 
calculations are based on the operator p s (x) = p 2 (x) (and smeared versions thereof). The 
corresponding (normalised) susceptibility is 

V((PS 2 } 1 



x aH</>.) a 

where V is the lattice volume. Fig. g(a) shows how this susceptibility varies with k at 
(3 = 6.0 and A = 0.01. While there is a pronounced peak, it neither grows nor sharpens 
as V is increased (except perhaps for very small volumes). Fig. 0(b) shows a similar plot, 
with similar behaviour, taken at (3 — 9.0. In Fig. 0(c) we show how the lightest scalar mass 



varies across this peak. (This is calculated from smeared versions of p s ; the smearing is 
exactly as in except that fundamental parallel transporters are replaced by their adjoint 
homologues.) According to a general finite size scaling analysis |7j a true second-order 
transition would be signalled by a peak in the susceptibility growing as Xmax ~ V 1 with 
some critical exponent 7. Likewise, across a second-order transition the lightest mass should 
go to zero with increasing lattice size, corresponding to a diverging correlation length. In 
the vicinity of the peak of the susceptibility we found the lightest scalar to be the lightest 
excitation. Although we do observe a susceptibility peak along with a minimum in the 
scalar mass, neither shows the required finite-size behaviour for a second-order transition. 
Instead we find \i «M[0 + ] — > const for V — > 00, which is the expected behaviour for 
crossover. The same observation holds at still larger values of A and at (3 = 9.0 as well. 

Where the first-order transition ends and the crossover begins, one would expect to find 
critical behavior, i.e. diverging fluctuations and a vanishing mass for the lightest scalar. 
This would form a line of critical points in our three dimensional coupling constant space. 
To investigate this, we calculate the minimum value of the mass of the scalar as we traverse 
the crossover, and then move along the line of minima towards the first-order region. 
Suppose we work at fixed A and j3, and find the minimum scalar mass m min (f3,X) at 
K — K min{P, A). If we now steadily decrease A at fixed (3, a critical point A c (/3), may occur 
at the boundary between crossover and first-order regions. This would be signalled by a 
vanishing of the scalar mass: m min oc (X — X c y. We have performed such a study at (3 — 6.0 
in which we have calculated m min for three values of A. Since m min is small, it is crucial 
that we ensure that the lattices we use are large enough. Our scans through the crossover 
were carried out on a range of lattice sizes from 8 3 to 32 3 . 

The results of these calculations are shown in Table |l|. We observe that m min (X) does 
decrease by about a factor of 2 when A is reduced from 0.011 to 0.009. This is suggestive 
of a critical point at a value of A < 0.009. Recall that we have performed calculations at 
j3 = 6.0 and A = 0.007 which reveal a very weak first-order transition. This would place 
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A 




d7n m i n 


0.011 


0.39820 (30) 


0.207 (23) 


0.010 


0.39600 (25) 


0.161 (14) 


0.009 


0.39380 ( 5) 


0.116 ( 8) 



Table 1: Minimal scalar mass in the crossover region approaching A c at f3 — 6.0. 

the suspected critical point for (3 = 6.0 in the interval 0.007 < A c (6.0) < 0.009. To locate 
it with more precision would require long runs on large lattices. 

We consider now gauge bosons and magnetic monopoles. Where the scalar particles are 
in the fundamental representation of SU(2), the existence of a residual or custodial SU{2) 
global symmetry allows operators which couple directly to the triplet of gauge bosons to be 
constructed and the gauge boson masses determined from correlators of these operators Q. 
For the present case of adjoint scalars, no such residual symmetry exists, and consequently 
there are no gauge-invariant operators which couple directly to charged gauge bosons. A 
gauge-invariant operator which directly couples to the photon in the Higgs phase can be 
constructed, however, starting with the continuum operator 2Tr(<pF flu ). In the Higgs phase, 
A 3 ^ is the electromagnetic gauge field in unitary gauge (0 = 3 o"3/2). In this gauge our 
operator is ~ F 3 U , and the latter is ~ d^A 3 — d u A 3 l , since the gauge fields A 1 ^ 2 have only 
short range fluctuations. Thus in the Higgs phase, 2Tr(<f)F f j -u ) couples to the field strength 
of the photon, and we can determine the photon mass by measuring its correlator ||. 
Working as we do in momentum space, at least one of the momentum components p M , p u 
must be non-vanishing in order to retain that part of 2Ti(<f)F^ u ) which is ~ d^A 3 — d^A^. 

On the lattice, we start from J 7 ^ = 2Ti(§(x)U flu (x)) which contains 2Tr(<pF tlu ) in the 
continuum limit. The operator we use in practice has the field $ inserted at all corners of 
the plaquette to maintain symmetry, and is then smeared to improve overlaps. We have 
measured the energy of the lowest state coupled to this operator at finite momentum. 
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Figure 3: (a) Photon energy for p = 2n/L at (3 = 6.0, A = 0.01. The solid lines mark 
p = 2tt/L for the respective lattice sizes, (b) Dispersion relation for the photon at (3 = 
6.0, A = 0.01, k = 0.4. The solid line marks the expectation for a massless photon. 

Fig. Q(a) shows the energy, E(p), for the lowest non-zero lattice momentum, p^ = 
17(1, 0) taken through the crossover at j3 = 6.0, A = 0.01. As can be seen, once the value of 
k is in the Higgs phase, i.e. k > k c (Q.0, 0.01) ~ 0.396, the energies agree very well with 2tx / L 
which is the energy of a massless excitation on a lattice of side L. The results are equally 
good for all of the three lowest non-zero lattice momenta that we studied. To demonstrate 
this we plot in Fig. |](b) E 2 (p) versus p 2 for all L and for the three lattice momenta at each 
L at a point well inside the Higgs phase. We see very striking evidence for the relativistic 
dispersion relation of a massless photon; a fit to the data yields (aM 7 ) 2 = 0.0010(20). In 
contrast, as we decrease k and move into the symmetric phase, the energy rises rapidly as 
it should. In the symmetric phase the operator in question describes an object composed 
of adjoint scalars and gauge particles or 'gluons' which will have a mass on the hadronic 
scale. Equivalent results are obtained at (3 = 9.0 and other values of A. 
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A photon of negligible mass in the Higgs phase is at first sight surprising. The classic 

work of Polyakov || established for adjoint scalars coupled to SU(2) gauge fields, that the 

photon has a finite mass due to the presence of 't Hooft-Polyakov monopoles. This apparent 

paradox is resolved by noting that for the present range of parameters, monopoles have a 

very large action and the photon mass established by Polyakov is far too small to show up on 

the lattices studied here. An estimate for the action of a monopole is S m ~ AnMw/g 2 , where 

Mw is the mass of the charged vector boson in the Higgs phase ||. We have not directly 

measured M w as this would require gauge fixing (see the discussion above). We proceed 

by using the tree level expression M^/ g A = —^rn^j (fig 2 ). Converting our parameter values 

(3 = 6.0, A = 0.01, k ~ 0.4 to the continuum using Eq.(|) we find S m ~ 5.8. Now Polyakov's 

formulas for the photon mass M 7 and monopole density n are M, 

w 2 2 {M w )l 

M -n, n exp(-5 m ). 

1 9 g 

Evaluating these at (3 — 6.0 gives (M 7 a) 2 ~ 10~ 4 and na 3 ~ 10~ 4 , so our small value for 
the photon mass is to be expected. We have measured the monopole density (in unitary 
gauge) using standard methods |10| in our configurations and find that there are indeed 
very few monopoles in the Higgs phase, as is displayed in Fig. |](a). Of those that are 
there, all (within errors) are members of tightly bound dipoles, and therefore ineffective in 
generating a photon mass. The signal for dipoles is that the average flux from a monopole 
vanishes into noise within one or two lattice spacings away from the monopole location. 

If, for all practical purposes, the Higgs phase does not contain any magnetic monopoles 
then as we cross over from the symmetric phase to the Higgs phase, the string tension 
(the coefficient of the linear piece of the confining potential) should vanish: a — > 0. To 
investigate this we have calculated the mass of the lightest flux loop that joins onto itself 
through the periodic boundary of the lattice. If we have linear confinement then the mass, 
mi, °f such a fl ux 1°°P should satisfy 

ami(L) = a 2 aL - — , (6) 
6L 
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Figure 4: Properties of the crossover region at (3 = 9.0, A = 0.01. (a) Number of monopoles 
in a 12 3 cube, (b) String tension; the upper solid line marks that in the pure gauge theory. 

where we have included the well-known universal leading correction to the linear piece. 
We calculate ami(L) from the asymptotic exponential decay of correlations of smeared 
Polyakov loops on L 3 lattices exactly as in the pure gauge theory. In addition to the 
standard Polyakov loops, lp, which we may schematically define by lp = Tr{]T^i U t }, 
we also consider modified Polyakov loops obtained by inserting scalar fields; i.e. lp = 
Tr{$(L)n t L =i^} and Zff = Tr{n t L =i U t $(t)}. These operators have the same properties 
under the centre of the gauge group as the Polyakov loop and therefore project only onto 
states with non-trivial winding, such as our periodic flux loop. We mention this technical 
detail because, as we move through the crossover region from the symmetric to the Higgs 
phase, the best operators become lp and then lp and one needs to use them rather than l P 
in order to obtain accurate mass estimates. (Of course, this observation tells us something 
interesting about the dynamics of the crossover as well.) 

We show in Fig. f|(b) a plot of our calculated value of the string tension, a 2 ex, as we vary 
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k through the crossover at /3 — 9.0 and A = 0.01. The calculation was performed as follows. 
We calculated mi(L) for L = 16,20,24 and, for perhaps half the K-values, for L = 32 as 
well. We then fitted these masses to the form Eq. (|6]). We obtained acceptable values of 
X 2 for k < 0.3805 and so for this range of k this is how we obtained a. (It is interesting 
to note that, throughout this range of k, fits without the 0(1/ L) string correction are 
very much poorer, even in the transition region where a 2 a is rapidly decreasing. This 
provides strong statistical evidence for the presence of such a correction and, indeed, for 
its 'universality'.) For k > 0.3810, the value of x 2 increased sharply. The confidence level 
of the best fit dropped from 90% at k = 0.3805 to 0.5% at k = 0.3810 and to < 0.01% at 
k > 0.3815. In this range of k we therefore allowed a constant piece and fitted the loop 
masses to ami(L) = c + a 2 aL. Such fits had an acceptable x 2 ■ (One can include a string 
correction as above, but it makes no significant difference.) 

We see in Fig. |](b) that the string tension collapses to zero around k ~ 0.3805. This 
is precisely the centre of the crossover (i.e. the coupling where the scalar mass attains its 
minimum value as we traverse the crossover at fixed (3 and A). For k > 0.3810, mi(L) 
is constant to a first approximation (apparently with a correction that seems to lead to a 
slight decrease with increasing L so that the fitted a has a tendency to be slightly negative). 
This is in fact what we would expect in a Higgs phase with no monopoles. The flux loop 
consists of a unit flux traversing the L x L spatial plane in, say, the y-direction. In the 
Higgs phase we would expect there to be no flux 'tube' so that the flux simply spreads 
out uniformly in the x-direction. Moreover, since we are in the [/(l)-like Higgs phase, we 
neglect non-Abelian effects. Then the field strength is E ~ the energy density ~ -j2 and 
hence when we integrate this over the Lx L plane we get a total energy that is independent 
of L. This matches what we observe and suggests that what happens to the confining flux 
tube as we traverse the centre of the crossover is that its width goes from being 0(-^) to 
being O(j^). This change takes place over a remarkably narrow range of couplings: between 
k = 0.379 and k = 0.381 in the case being considered here. Of course if we let L — > oo 
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then there will be a very dilute gas of monopoles and this will reintroduce a finite, albeit 
extremely wide, flux tube and a non-zero, albeit extremely small, a. Finally we remark 
that although the Higgs phase may have no visible string tension, it is still confining, albeit 
only logarithmically, because we are in two spatial dimensions. Thus (lp) = for all values 
of the coupling and there is no need for a phase transition. 

We return briefly to the lines of constant physics. At (3 = 9.0, k = 0.3737, A = 0.01 
we are deep in the symmetric phase. Using a variety of operators we measure masses and 
the string tension as above. Now along a line of constant physics the ratio of a mass to 
the square root of the string tension should be independent of the lattice spacing scale, if 
the 0(a) corrections are small. To test this, we use Eqs.([|) and (|5|) to map k = 0.3737, 
A = 0.01 to the corresponding values at (3 = 6.0, and see whether the two sets of calculated 
mass ratios agree. This is indeed what we find for operators of purely gluonic, purely scalar 
or mixed origin (within errors which are < 10%). Repeating this for a (3 — 9.0 point in the 
crossover region (A = 0.01, k = 0.3805), agreement is still good, although less so for the 
admittedly more noisy gluonic operators. Far inside the Higgs phase (A = 0.01, k = 0.3900) 
the string tension is vanishingly small. Mass ratios constructed using the lightest scalar 
mass (obtained from the purely scalar operator) as a reference do give good agreement for 
the operators of mixed origin, however. Thus we have good evidence for scaling between 
(3 = 6.0 and 9.0, along lines of constant physics, and hence that our conclusions in this 
paper are relevant to the continuum theory. 

In conclusion, we have clarified the phase diagram of the SU(2) adjoint Higgs model 
in 2+1 dimensions. For small values of the scalar self-coupling, A, the expected first- 
order transition indeed occurs. But as we increase A this becomes a crossover rather than 
a second-order transition. At the point where the first-order transition vanishes there 
appears to be a critical point where the scalar mass vanishes. These calculations have been 
performed for values of the lattice spacing where deviations from the continuum limit are 
very small and so we believe that these are properties of the continuum theory. All this is, 
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of course, very much reminiscent of the behaviour of the theory with fundamental Higgs 
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Understanding the phase diagram is an essential first step if we wish to use this theory 
as a testing ground for the monopole mechanism for confinement. The fact that the Higgs 
phase and the symmetric, confining phase are smoothly connected through a crossover 
region, makes us optimistic about the possibility of learning how confinement interpolates 
between the former phase, where it is well understood, ||, and the latter, where theoretical 
ideas exist, ||11|| , but still remain speculative. This optimism is reinforced by the simple 



behaviour that we have found for the flux tube and the string tension as we pass through 
this crossover region. Of particular interest in this context is the critical point; it may be 
that we can find here a reasonably dense plasma of t'Hooft-Polyakov monopoles. These 
questions are under investigation. 
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